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Abstract 

The Operator Product Expansion, in conjnnction with the power connting of non- 
relativistic held theory, is used to examine the end-point region of the radiative 
decay of heavy quark bound states with ^Si qnantum numbers, QQ 7 + 
X. We identify an inhnite class of operators that determine the shape of the 
photonic end-point spectrnm. These operators can be resummed to form an 
octet strnctnre fnnction which parameterizes the energy of the dynamical glnon 
content of the leading octet Fock state component of the quarkoninm. This color- 
octet contribntion is important when the photon spectrnm is examined with a 
resolution given by AE^ ~ mqv'^, where is the relative quark velocity squared. 
The formalism used makes explicit the shift of the end-point from its partonic to 
its hadronic value. 


The direct photon spectrum in T(IS') decay has been studied experimentally|]l[] and 
compared with predictions based upon the color-singlet model[^ |], ||, |]. These predictions 
were rehned by the inclusion of hadronization effects which one would expect to be 
important in the end-point region. In this letter we discuss an effect which is formally larger 
near the end-point which has been heretofore neglected. This effect stems from higher Fock 
states in the heavy quark bound state. Using an operator product expansion (OPE) in 
conjunction with non-relativistic QCD (NRQCD) |^, we show that the end-point spectrum 
receives a leading order contribution from a color octet structure function. Furthermore, the 
formalism shows explicitly the shifting of the maximal photon energy from its partonic to 
its hadronic value. Thus, studying the endpoint region of the spectrum becomes a useful 
laboratory for studying non-perturbative QCD effects. 

Direct photons in quarkonia decay arise from the electromagnetic coupling of both the 
heavy and light quarks. The spectrum stemming from the couplings to the light quarks 
has been previously studied in ^ using the leading logarithmic approximation. This frag¬ 
mentation contribution to the photon energy spectrum is suppressed relative to that coming 
from the coupling to heavy quarks, especially near the photonic endpoint where the hadronic 
invariant mass is highly restricted. In this letter, we will focus on the spectrum due to the 
coupling to the heavy quarks. 

The photonic decay of quarkonia has the merit that the decay rate contains an external 
variable, namely the photon energy, which allows one to analytically continue to the unphys¬ 
ical regime where the OPE is trustworthy. This is tantamount to the statement that we may 
smear over a range of hadronic invariant mass, resulting in an average which is dual to the 
parton model 0. This is in contrast to the total hadronic decay calculations which assume 
duality on a point by point basis. 

For the case at hand, we begin by considering the imaginary part of the time ordered 
product of two electromagnetic currents 

T = iJ I T(J^(a;) J^(0)) | T), (1) 

where = 67^6 and the normalization convention for upsilon states is (T(p) | T(p')) = 
(27r)^h^(p^—p). In eq.(l) and hereafter, spin averaging of the T matrix element is understood. 
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At fixed, = 0, T has two cuts which span the entire real go axis (we work in the rest frame of 
the T) in the complex go plane. There is a cut resulting from the physical process T —> 7 + X 
in the region —mx /2 < go < mT /2 (neglecting the pion mass) and another from the process 
7 + T —> X for which the hnal state, X, has invariant mass larger than mx- This overlap 
obstructs the usual dispersion relation analysis since the discontinuity across the cut in the 
physical region is polluted by the other process. However, this obstruction is a red herring. 
To see this, we note that the T —7 + X decay rate can be extracted from the time ordered 
product 

T' = tl I T(J^(x)i2 J^( 0 ) 34 ) | T34), ( 2 ) 

where J/)- = bi'y^bj and bi are differing degenerate b-type quark species. T' has two cuts, 
however they are now separated as shown in Figure 1. To get the T ^ 7 + X rate we simply 
perform the contour integral over a contour which only picks up a contribution from the part 
we are interested in. We may then relate the decay rate to the imaginary part of the T' via 


dr 

dE^ 


TT 


( 3 ) 


where = —1/3 is the electric charge of the 6 -quark and a is the hne structure constant. 

We now may perform an OPE on the product of currents and calculate the Wilson 
coefficients using perturbative QCD. The OPE is valid in the unphysical region which is 
related to the physical process via the contour shown in Figure 1. The contour approaches 
the cut at one point, and thus the level of rigor here is not at the same level as in deep 
inelastic scattering but should be considered at the same level as the predictions for semi- 
leptonic heavy meson decay p, ^]. Performing the contour integral over the photon energy 
corresponds to the smearing mentioned above. We may choose different weightings for the 
contour to extract information regarding the spectrum. Furthermore, as will be seen below, 
the convergence of the OPE will be dictated by the size of the smearing region. 

For a given photon energy E^, the hnal hadronic invariant mass is given by — 

2E^/mr)- Consequently, the end-point region of the photon energy spectrum corresponds 
to low invariant mass hadronic hnal states. In this region, the photon energy spectrum 
eq.(3) must be smeared over a range AE^. We will see that the OPE gives us quantitative 
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Figure 1: Analytic structure of T' in the complex go plane and the appropriate choice of 
contour. The cut on the left hand side begins at go = fn-x/2 whereas the cut on the right 
hand side begins at go = 3mT/2. 

information regarding the necessary size of the smearing region. Indeed, the convergence of 
the OPE is predicated upon the prudent choice of smearing functions. 

The leading order diagrams come from those shown in Figure 2b. This decay mode is 
dominated by the color-singlet Fock state. NRQCD velocity scaling rules place this contribu¬ 
tion to the decay rate of order, Away from the end-point this is the leading piece. 

However, near the end point there will also be a contribution stemming from those diagrams 
shown in Figure 2a which are dominated by the Fock state where the quarks are in a relative 
octet state. This color-octet process yields contributions to the decay rate of order aa^v^ jn 
and aagv’^/n for the ^S'o and ^Pj intermediate states, respectively [^. If we use ~ Os/vr, 
then the contributions are down by a factors of v and relative to the color-singlet process. 
However, this naive counting does not hold for the end-point region of the photon energy 
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spectrum because the color octet contribution is highly singular there. Indeed, as will be 
shown below, after smearing over a region of photon energies, AEj ~ rubv"^, the octet process 
is of the same order as the singlet in the endpoint region. 

The leading order octet contribution to T' is given by 




+ i + mb) 
{p + iD — qY — ml 


^T^lub2 


{2p — q + iVy + ie 


Tj,A + + 

^3-^ 1 U / I ' 7 ^ NO 0 4 . 

[p + %D — qy — mi 


T 34 ), 


where each of the heavy quarks carries momentum p = {mb, 0) and D denotes a covariant 
derivative. Tg receives contributions from ^5*0 and ^Pj intermediate states so we write 


T' = T'{^Pj)+T'CSo). 


We now match onto the NRQCD by expanding the helds b in terms of non-relativistic 
helds Ip and y for the particle and anti-particle respectively. In making the transition to 
the effective theory, we keep only the Dirac strnctnres relevant to the decay at hand. The 
relation between b-qnark helds in fnll QCD and NRQCD prodnces additional factors of D. 
After imposing spin symmetry we hnd that the contribntion to the imaginary part of T' 
from the octets Fock state are given by 

/mT'CSo) = glCis, J dk*OCS„)f{kp,sAmi. - A + (6) 


IrnTAPj) = glChp, I dk+0{Aj)f(kp,pAmt - E, + k* 12), (7) 

where we have dehned 

D, TA\ I T)IO('Pj), (8) 

0(Aj) = (T I Dj [xV.j D, Ap I T>, 

fiKhs, = (T I S(k+ -n-W) [x'tA] I T)/0(‘S„), 

0(‘S„) = (T I [x'T-V] I T), 


flKhpj = A \ S(k^- n-iD) x*fip 
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(9) 

( 10 ) 

( 11 ) 



and the light-like four-vector n = {2p — q)/mb is taken to be, n = (1, 0, 0,1). The four-vector, 
kg , is the light cone momentum of the dynamical gluon in the octet Fock state. Furthermore 
we have kept only the leading twist operators, which is to say we have dropped terms of 
order {kg)^ in the delta function. The constants Capj and have the values 


C^pj 


Ttt 

18mf ’ 


a 


So 


71 

2 mf 


( 12 ) 


The functions f{kg)3pj and f{kg)isQ are the normalized probabilities that the dynamical 
gluon in the octet Fock state has light-cone momentum kg . A strong analogy may be made 
with the deep-inelastic structure functions by studying the Fourier transform of the /(/cj)’s. 
In the gauge, n • A = 0, the Fourier transform introduces a factor of exp{tn ■ d). This 
translates the combination of quark helds [x^'0] from the origin to the spacetime point tn 
giving, for example in the case of ^Pj, 


= (T 


Dj T^x(O) 


Pexp[—igs fo dt'n ■ A{t'n)] 

om) 


X^{tn)ai- Dj 




T). 


(13) 


In eq.(13) a path-ordered exponential has been inserted to restore gauge invaiance. 

The analysis of the color singlet contribution is simplihed by the fact that in this case 
we may use the vacuum saturation approximation, which is a well controlled expansion in 
for this case 0]. The net result of summing the leading twist corrections is simply to shift 
the maximal photon energy from its partonic to hadronic value. To see this, we note that 
the imaginary part of Figure 2a can be written as 


ImT^Si) = j;G(£),,) y 


~ I | T), 

^ n! d[2E^D ' L j L j 

(14) 

where the leading order Wilson coefficient, G{E^), is a smooth function of E^, that does 
not vanish at the end-point. After factorizing this matrix element only the time derivatives 
contribute (for T states at rest) giving factors of the binding energy |]^. Consequently 
eq.(14) becomes 


ImTiCSi) = gtG{E,)e{mr - 2P,)(T | I 0)(0 | x^a.i^ \ T), (15) 
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2b. 

Figure 2: 2a) Leading order diagram for the octet contribution. 2b) Leading order Diagram 
for the singlet contribution. 


showing that the whole effect of the leading twist contributions for the color singlet part of 
the amplitude is to shift the end-point from its partonic value, nib, to the physical one m-r/2. 

We are now in a position to consider the relative sizes of T/(^S'i), T^i^Pj) and Tg(^S'o) 
in the end-point region. If the smearing region is large enough, then we would expect that 
the singular nature of the octet contribution will not compensate for naive v suppression of 
the octet state contributions. Too small of a smearing region will lead to a breakdown of 
the OPE, as the subleading twist terms will begin to dominate the octet contribution. More 
quantitatively, we may expand the octet contribution into an inhnite sum of leading twist 
matrix elements 


ImTi(nj) = C.p.gJ 

n=0 



Dj T\ 





T) (16) 
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/mT'(‘S„) = C\s„3l Y. W-{T I V'T-^x 




n=0 


X^T^p] I T) (17) 


where <5^"^ denotes n derivatives with respect to 2E^ acting on the delta function. Eqs.(16,17) 
are similar to the expression found in the analysis of the leptonic end-point spectrum in B 
decays P!T| , |^. An operator with n factors of the covariant derivative scales as [|^ y(6,7)+2n^ 
as each derivative contributes a factor of to the scaling. If we smear, for example, using a 
Gaussian of width AE^ then for, AE^ << the terms in the sum grow with n, and the 

subleading-twist terms which have been dropped cannot be neglected. On the other hand, for 
AEy » the singlet contribution will dominate over the octet. Thus, the hnest energy 

resolution with which we can examine the photon endpoint spectrum, without introducing 
yet more higher twist structure functions, is of order rUhV^. Note that this corresponds to 
smearing over a range of hadronic masses Amx ~ "^bV which for large is much greater 
than the QCD scale. For the smearing width AE^ ~ rubv"^ the octet ^Pj contribution is of 
the same order as the color-singlet since the naive octet suppression of is compensated by 
the fact that the singlet contribution starts with a theta function, whereas the octet starts 
with a delta function. The formally leading contribution Q however, comes from the octet 
^S'o contribution which is enhanced by a power of 1/v. Therefore, if we wish to resolve the 
photon energy spectrum within rribv'^ ~ 500 MeV of the end-point, we must know the value 
of the structure functions f{k'^)isQ and f{kg)3pj, which can be perhaps measured by the 
lattice or extracted from another process. 

Let us now consider the size of other possible large effects near the end-point. As men¬ 
tioned in the introduction there is another nonperturbative effect near the end point having 
to do with the hadronization of the partons. The gluonic partons are not really massless, 
but are actually hnite invariant mass objects. This effect was addressed in ref. p[, where 
the author used a parton-shower Monte Carlo approximation to determine the effects of 
hadronization on the photon spectrum of the color-singlet contribution. Hadronization cuts 
off the photon distribution at the end point. This effect is important for hadronic masses. 


^ See erratum in ref. |^. 
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mx, of order the QCD scale and should be subdominant to those of the octet contribu¬ 
tion. The success of hadronization models in describing the endpoint spectrum may indicate 
that the color octet matrix elements are smaller than expected on the basis of dimensional 
analysis. 

In the end-point region there are also large perturbative corrections stemming from terms 
of the form ln(l — x) where x = E^/rrib- These terms result from the fact that near the end¬ 
point gluon radiation is suppressed ruining the delicate cancelation between real and virtual 
graphs. A resummation of these logarithms is expected to suppress the singlet spectrum 
near the end point. The calculations in ref. [Q indicate that this resummation does indeed 
suppress the color-singlet rate near the end-point and are imporant when > O.dmx- 
There will be similar large logarithms in the Wilson coefficient for the octet contribution 
discussed above, resumming these logarithms should again suppress the end point region. A 
calculation of these effects is necessary to undertake any phenomenological investigation of 
the end-point region. 

The results of this paper are applicable to both cc and bb quarkonia states. However, 
we feel that the charm quark mass is probably not large enough for our methods to be 
trustworthy. We also note that effects similar to those discussed here are important for 
fragmentation into quarkonia near z = 1 a.s well as direct production in some kinematic 
circumstances [jl^ . 
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